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NONLINEAR OBSERVER IN THE ORIGINAL COORDINATES WITH 
DIFFEOMORPHISM EXTENSION AND JACOBIAN COMPLETION 

PAULINE BERNARD, VINCENT ANDRIEU AND LAURENT PRALY * 


Abstract. Typical difficulties we face in implementing observers in applications are: 

- because of physical constraints or observability singularities, the estimated state should remain in some restricted 
region of the system state manifold ; 

- observability may also impose the observer dynamics to live in an observer state manifold with higher dimension 
than the system state. In such a case, the observer implementation needs a left inverse of an injective immersion. 

In the approach we propose to round these difficulties, we do not modify the given converging dynamics of the 

observer, assumed to be globally defined. Instead we select carefully the coordinates to express them. This is done 

via appropriate diffeomorphisms: 

- one diffeomorphism sending the restricted region to the whole space. This is done by an image extension. 

- one function extending the injective immersion into a diffeomorphism and making the observer state coordinates 
equivalent to the system state coordinates complemented with extra ones introduced to fill the dimension gap. 
This is done via a Jacobian completion. 

Such a design makes possible the expression of the observer dynamics in the, maybe complemented, original coordi¬ 
nates. 

Several examples illustrate our results. 


1. Introduction. 

1.1. Context. In many applications, estimating the state of a dynamical system is crucial 
either to build a controller or simply to obtain real time information on the system. A lot of efforts 
have thus been made in the scientific community to find universal methods for the construction 
of observers. Although very satisfactory solutions are known for linear systems (ini), nonlinear 
observer designs still suffer from a significant lack of generality. 

For nonlinear systems, we are aware of two ’’general purpose” observer design methodologies 
with guaranteed ’’non local” convergence: the high gain observers ( [161 [2llHI[9] etc) and the non¬ 
linear Luenberger observers ([201 Ha El). Those observers do not demand any particular structure 
and only require some basic observability properties. However, in both cases, the observer state is 
in general living in a space different from the system state one and the state estimate is obtained 
typically by solving on-line a nonlinear equation, which may be very complicated. 

As an illustration, consider an harmonic oscillator with unknown frequency with dynamics 


±1 = X2 , ±2 = -X1X3 , ±3 = 0, y = Xi 


( 1 . 1 ) 


with state x = (xi, X 2 , X 3 ) in x IR>o and measurement y. We are interested in estimating as x 
the state x from the only knowledge of y and maybe the fact that x evolves in some known set. We 
can solve this observer problem by following in a very orthodox way the high gain observer design 
(see [15] for example for the general theory and [1] for details). This leads to an answer with a 
dynamical system, with dynamics 
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with observer state ^ in and where sat is a saturation function. With this, the system state 
estimate x is obtained as 


X = T{i) 


where the function r is obtained from solving in x 

1 = (ll, 12,13,14) = T*(x) = (xi,X2,-XiX3,-X2X3) 


(1.3) 


This is a system of 4 equations in 3 unknown which in general has no exact solution. To get an 
approximate solution, we may have to solve an optimization problem such as 


X = t{^) = Argmjn ^ — t*{x) 


We are aware that, for the particular expression of t* in (1.31, an expression for t(^), solution of 
the above problem, can be obtainecj^ But this is not the case in general and we go on in this paper 
ignoring the data of r to get x from Actually we cannot fully ignore r since the values of (p in 
(1.2) depend on x. This subtlety will be addressed in Assumption A.C below with dealing with 


pairs {ip, t). 

To eliminate the minimization step above, we may want to write the observer dynamics directly 
in the x coordinates. This is not straightforward since x has dimension 3 whereas ^ has dimension 
4. We round this difficulty by adding one component, say w to x. With this, we are left with 
writing the dynamics of (x,w). This can be done by moving, in the {x,w) coordinates, the f 
dynamics known to solve the observer problem. For this, we nee d a diffeomorphism from [x, w) 
to ^ which should “extend” the function x 1 —)■ t*{x) given in (1.3) We show in Section that this 
can be done by a Jacobian completion. Unfortunately, in doing so, the obtained diffeomorphism is 
rarely defined everywhere. Thus, we are facing the new problem of having no guarantee that the 
trajectory in {x,w) of the observer remains in the domain of definition of the diffeomorphism. We 
show in Section]^ how this difficulty can be rounded via a diffeomorphism extension. 

In the following, as in the example above, we assume we are given a preliminary, say raw, 
observer with state ^ of dimension possibly different from the one of x and we look for a more 
satisfactory one written in the, may be complemented, x components. 

Writing the observer in the original coordinates has been suggested by several researchers 
0 na la in the case where the dimension of the raw observer state and the system state are 
the same. Our contribution is to relax this latter condition via Jacobian completion for which 
preliminary results have already been presented in [1]. 

As mentioned above in the example, the route via Jacobian completion may suffer problems 
with the domain of validity of the diffeomorphism. Actually, we may have such a validity problem 
already with the raw observer (see Section]^. This makes necessary the study of diffeomorphism 
extension which in the context of observer designs has not been studied yet, as far as we know. Its 
objective is to force the observer state to remain in a specific set without modifying the given raw 
observer dynamics. Another path, which has been proposed in UHlEl, is to modify the dynamics. 
But then, extra assumptions such as convexity are needed to preserve the convergence property. 


^Another possible expression is t(|) 
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We motivate and illustrate our results with continuing the example of the harmonic oscillator 
with unknown frequency and with the bioreactor presented in [5]. These are done with using a 
high-gain observer as raw observer only, due to space limitations. But, as shown in [4], exactly the 
same tools apply to nonlinear Luenberger observers a la uni in HI- 

1.2. Problem statement. We consider the given system with dynamics : 

i = f{x) , y = h{x) , (1.4) 

with X in K" and y in K'^. The observation problem is to construct an algorithm which, from the 
knowledge of y, estimates the system state x as long as it is in a specific set of interest denoted 
A C R”. In this paper, our starting point is to assume that this problem is (formally) already 
solved but with maybe some implementation issues. More precisely, we assume : 

A.II=Injective immersion : There exists an open subset O of M" which contains A and a 
injective immersion t* : O ^ R™. 

Actually this injective immersion should be such that the image LfX* by r* of the vector field / , 
is in a form such that an observer can be designed. We write the corresponding observer as : 

i = (p{i,x,y) , x = T{i), ( 1 . 5 ) 


where the functions (p : R™ x R" x R'J —> R" and r : R™ —> R" are locally Lipschitz. For the image 
by T* of a solution of (1.41 to be a solution of (1.51, at least as long as the former is in A, we must 
have : 


LfT*{x) = ip{T*{x),x,h{x)) 


^{t*{x)) = X Vx e A 


( 1 . 6 ) 


In the following we do not impose explicitly these identities are satisfied. They are imposed implic¬ 
itly via the convergence requirement made precise below. 

As already noticed with (1.2), we need the function r to reconstruct x, argument of in (1.51. 
This dependence of ip on r may imply to change p whenever we change r. This leads us to consider 
a se10</^ of pairs {p,t). With this, we can now phrase our assumption concerning the convergence 
property mentioned above. 

A.C=Convergence : We know a set AT of locally Lipschitz functions {p,t) such that, for any solution 
X{x, t) of { 1 . 4 ) which is defined and remains in A for t in [0, -foo), the solution {X{x, t), S((^, x),t)) 
of the cascade system : 


issued from (xA) in Ax 


i = fix) , i = pii,T{i),hix)) , 

at time 0, is also defined on [0, -foo) and satisfies : 


(1.7) 


lim 

)- + oo 


T*iXix,t)) - ^iii,x),t) 


= 0 . 


Very trivially, p can be paired with any function r in the particular case in which p does not 
depend on r, i.e. it is in the form pi^, t{x), y) = p{^, y) as in the case of the nonlinear Luenberger 


^The symbol i^is pronounced phitau. 
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observer for instance (see [1]). In the high-gain approach, when A is bounded, thanks to the gain £ 
which can be chosen arbitrarily large, (p can be paired with any locally Lipschitz function t provided 
its values are saturated whenever they are used as arguments of (p. 

In the following, we consider that the function r* and the set (pT are given and we aim at 
designing an observer for x in the, maybe complemented, x-coordinates. However, in trying to 
meet this objective, we shall encounter some difficulties which will lead us at the end to reconsider 
the data of t* and (/ST 


Example 1. For the harmonic oscillator with unknown frequency (1.1), since, for any solution 


with initial condition xi = X 2 = 0 , we cannot get any information on X 3 from the only knowledge 
of its dynamics and the function t 1 —>■ y(t) = we define the set A in which we want the 

estimation as 


A =< X G 


xj 


xl G 


1 

-,r 

r 


, X 3 G]0,r[ 


( 1 . 8 ) 


where r is some arbitrary strictly positive real number. This set is invariant and, the system is 
stron gly differentially observable of order 4 on this set. Indeed, it can be checked that the function 
(1.3) is an injective immersion on (M^ x M>o) \ ({(0,0)} x K>o). This implies that Assumption 
A.II is satisfied for any open subset O such that cl(A) G O C (K^ x K>o) \ ({(0, 0)} x K>o), with 
cl denoting the set closure. Then, relying on what is known on high gain observers (see |15j for 
example), we can claim that, in the present case where A is compact, a set ipT satisfying Assumption 
A.C is made of pairs of a locally Lipschitz function r satisfying 


X = r(xi,X 2 , —X 1 X 3 , —X 2 ,X 3 ) Vx G A 

and the function (p defined in ( | 1 . 2 [ ) where 

sat(s) = min {r^, max {s, —r^}} (1-9) 

and the gain is adjusted to the properties of t. 

We do not need to know an expression for r, in what follows knowing its existence is sufficient. A 

1.3. A sufficient condition allowing us to express the observer in the original coordi¬ 
nates. To motivate the technicalities we shall use to write the observer in the original coordinates, 
we start with the simpler case where the raw observer state ^ has the same dimension as the system 
state X, i.e. m = n. The example of the bioreactor from [ 8 ] detailed in Section falls into this 
category. In this case, t* is a diffeomorphism on O, i-e it is simply a change of coordinates. We 
can implement the observer in the original coordinates : 

^=(^(^)) V>{'r*{x),x,y) ■ ( 1 - 10 ) 

which only relies on a Jacobian inversion. However, although we know by assumption that the 
system trajectories remain in O where the Jacobian is invertible, we have no guarantee the ones of 
the observer do. Therefore, we must find means to ensure the estimate x does not leave this set 
in order to obtain convergence and completeness of solutions. To achieve this, we will see that a 
possible approach is the extension of the image set t*{0) to K™. 

With this in mind, we now look at the more complicated situation where m > n, i-e t* is only an 
injective immersion. The example of the harmonic oscillator with unknown frequency corresponds 
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to this case. In [T], it is proposed to augment the given injective immersion r* into a diffeomorphism 
T* : Oe = O X Sw ^ R™, thus adding m — n dimensions to the state through a new variable w. To 
help us in finding such an appropriate augmentation, we have the following sufficient condition. 
Proposition 1.1. Assume A.II and A.C hold and A is bounded. Assume also 

PI: Completion to a diffeomorphism. There exists an open subset Oe o/K’" and a diffeomor¬ 
phism T* : Oe —)■ such that Oe contains A x {0} and we have 


T*{x,0) = T*{x) \/x € A . 


( 1 . 11 ) 


P2: The function 


is surjective. In other words, Tf{Oe) = R"*. 


PS: Couple {x,{p) admissible. With denoting Tex the x-component of Te, the inverse ofrf, there 
exists a function tp such that the pair {tsx, ff) is in the set ffTgiven by assumption A.C. 

Under these conditions, for any solution X{x,t) of which is defined and remains in A for 
t in [0,+oo) the solution (X{x,t), X{x,x,w,t),W{x,x,w,t)), with initial condition {x,w) in Oe, 
of the cascade of system (!./([) with the observer : 



drt 


d{x, w) 



-1 

T{Te{x,w),x,y) 


is also defined on [0, +oo) and satisfies : 


lim 

t —^-|-oo 


W{x, X, w, t) 


X(x, t) — X{x, X, w, t) 


= 0 . 


( 1 . 12 ) 


(1.13) 


The key point in the observer (1.121 is that, instead of left-inverting the function r* via r as 


in (1.5), we invert only a matrix. 

Proof. For any compact set £ containing A x {0} and contained in Oe, there exists a class JC 
function a : IR>o —t R>o such that we have 


\w\ \Xa-Xb\ < a{\T*{Xa) - Tf{xb,w)\) ^ Xa & A , ^ {Xb,w) & . 


(1.14) 


Indeed, let: 


ao{s) = max \xb - x^ + |w& - We\ 

{{xb,Wb), {Xc,Wc)) e 

\Tf{xb,Wb) - Tf{Xc,We)\ < S 

This defines properly a non decreasing function satisfying : 

\Xb - Xc\ + \Wb - We\ < ao{\Tf{Xb,Wb)-Tf{Xc,We)\\) V {{Xb, Wb), (Xc, We)) G ■ 

If q;o( 0) is not zero, then there exists {{xb,Wb), {xc,Wc)) in satisfying : 

\Tf{xb,Wb)-Tf{Xe,We)\ = 0 , \Xb - Xe\ + \Wb - We\ > 0. 


^ In the particular case where a subset S of is known to be invariant by the dynamics ^ r(j), y), it is 

enough to have S C T*{Oe) in P2. But dealing with this case is difficult when S depends on r. 
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But this contradicts the injectivity of t* (i.e. PI). Then let Sk be a sequence converging to 0. For 
each k and I in N>o, we can find {{xb,ki,Wb,ki)j ixc,ki,Wc,ki)) which satisfies 

^oi^k') ^ \^b,kl “t“ \'^b,kl 1 ^ ^o('^fc) y : 

I'Tg {Xb^kl ^'^b,kl) {Xc^kl ^ ■ 


The “diagonal” sequence {{xb^kkiWb^kk), (xc^kk^Wc^kk)) being in a compact set admits an accumu¬ 
lation point {xb,*,Wb,*), {xc,*,Wc,*) which, because of the continuity and the injectivity of r* must 
satisfy (x;,,*, Wb,*) = (xc,*, Wc,*)- This implies that ao(s/c) tends to 0 and thus ao is continuous at 
0. Now, consider the function defined by the following Riemann integral 

1 

a(s) = - ao{v)dv + s . 
s Js 

It is continuous, strictly increasing and zero at zero and we have : 


\xb- Xc\ + \Wb- Wc\ < a{\Tl{xb,Wb) - tI{Xc,Wc 


\/{{xb,Wb), {Xc,Wc)) e 


So in particular for Xc = Xa and Wc = 0, with Xa in A and (1.11), we obtain (1.14). 


Let (x,(x,w)) be arbitrary \n Ay. Of. but such that X{x,t) solution of (1.4) is defined and 
remains in A for t in [0,-l-oo). Let [0,r[ be the right maximal interval of definition of the so¬ 
lution {X{x,t),X{x,X:W,t),W{x,XtW,t)) when considered with values in x Oe- With aim¬ 
ing at showing that T is infinite by contradiction, assume it is finite. Then, when t goes to 
T, either (X(a:, ri), t), lF(a:, a;, w, t)) goes to infinity or t o the boundary of Oe- By construction 
(^X{x, w, t),W{x, w, t)^ is a solution of (1.7) on [0, T[. From P3 and assumption A.C 

= T:(Oe). 


considered with values in 


t E{t) := ■ 

it can be extended as a solution defined on [0, -|-oo[ when 
This implies that E{T) is well defined in K"*. Since, with P2, the inverse Te of r* is a diffeomorphism 

defined on M™, we obtain limi_>r (^X{x,w,t)^W{x,w,t)^ = T'e(S(F)), which is an interior point 

of Te)®™) = Oe- This point being neither a boundary point nor at infinity, we have a contradiction. 
It follows that T is infinite. 

Finally, with assumption A.C, we have : 


lim 

t —^ “|“00 


*{X{x,t)) — T* ^X{x,w,t),W(x,w,t)j 


= 0 


Since T*{X{x,t)) remains in the compact set t*{A), there exists a compact su bset C of M™ and a 
time tc such that r* u), t), lF(x, re, t)^ is in C for all t > tc- Applying (1.14) to £ = Xe{C) 

which is a compact subset of Oe, we obtain (1.13). □ 

Remark 1. It follows from this proof that the assumptions of boundedness of A and \1-11 ) 
can be replaced by: For any compact subset £ of Oe, there exists a class K. function a : ®>o —>■ K>o 
such that, (1.14) holds. Also conversely, if (1-14) o,nd P2 hold, then (TTTi holds. 

Addressing problems PI and P2 which have their own interest outside the observer context is 
the main topic of this paper. We will present in Section conditions under which the extension of 
an immersion into a diffeomorphism in the sense of PI can be done via a Jacobian completion. We 
will address problem P2 in Section with giving conditions under which we can extend the image 
of the diffeomorphism r* to M™. 
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Throughout Sections 2][3 we will show how, step by step, we can express a high gain observer in 
the original coordinates for the harmonic oscillator with unknown frequency. The various difficulties 
we shall encounter on this road will be discussed in Section]^ In particular, we shall see how they 
can be rounded thanks to a better choice of r* in Assumption A.II and the pair (r, Lp) given by 
Assumption A.C. Finally, we will show in Sectionj^that our approach enables to ensure completeness 
of solutions of the observer presented in [5] for the bioreactor. 


2. About PI : Jacobian completion. 

2.1. A way to satisfy PI. In [T], we have the following sufficient condition for having PI 
satisfied. 

Proposition 2.1 m )- Assume A is bounded and A.II holds. Let be a bounded open set 
such that A C cl{Sx) G O. If there exists a function 7 : cl{Sj;) —t ") whose values are 

m * (m — n) matrices satisfying : 


det 



j{x) 0 


Va; G cl{Sx) , 


( 2 . 1 ) 


then there exists a strictly positive re al n umber £<> such that, for any e in (0,eo], a function r* 
satisfying condition PI of Proposition [iT7| = Sx x Bs{0) is : 

Tf{x,w) = T*{x) + 'y{x)w . (2.2) 


Proof. The fact that r* is an immersion for e small enough is established in [T]. We now prove 
it is injective. Let £0 be a strictly positive real number such that the Jacobian of r* {x, w) in (2.2) is 
invertible for any {x,w) in cl(5a; x ^^^(O)). Since cl(5a; x Beg{0)) is compact, not to contradict the 
Implicit function Theorem, there exists a strictly positive real number S such that any two pairs 
{xa,Wa) and {xh,Wb) in cl(52; x ,Beg(0)) which satisfy 


*{Xa,Wa) = Tf{xb,Wb) , {Xa,Wa) ixb,Wb) 


(2.3) 


satisfies also 


\Xa - Xb\ + \Wa - Wb\ > S . 

On another hand, since r* is continuous and injective on cl(iSa;) C O, by following the same 


arguments as those for establishing (1.14), we can prove the existence of a class K, function /3 such 
that we have 

\Xa-Xb\ < P{\T*{Xa) - T*{xb)\) \/{Xa, Xb) € cl{Sx)‘^ . 


It follows that, if (2.3) holds with Wa and Wb in ^^(0) with £ < £ 0 , we have 


5 - 2e < |a;a-a;&| < l3{\T*{xa) - T*{xb)\) = l3{\^{xa)wa - ^{xb)wb\) < ,5 2£ sup | 7 (a;) 

\ a:Gcl(Sa:) 


■^For a positive real number e and zq in RP, is the open ball centered at zq and with radius e. 
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But there exists a strictly positive real number £<> < Eq such that, for any e in (0,eo]) we have 


S - 


2e > j3 


( 2 e sup 

xGcl(5a,) 



So, for all such e, (2.3) is impossible. 


This proves that r* is injective on Sx x Se( 0) for any e in 


(0,e«]. □ 

With this proposition, condition PI is related to the existence of a function 7 such that the 
matrix (13 is invertible. 


2.2. Submersion case. We can solve the Jacobian completion problem when T*(cl(5a:)) is 
a level set of a submersion. 

Proposition 2.2 (Completion from a submersion). Assume A.II holds. Let Sx be a hounded 
open set such that cl{Sx) C O. Assume there exists a function F : K™ —)■ R'"”" which is a 
submersion at least on a neighborhood of t*{Sx) satisfying: 


F{t*{x)) =0 \/x € Sx 


(2.4) 


then with the function x 1 —>■ 7 ( 3 :) = ^ {t*{x)), the matrix in (2.1) is invertible for all x in 


Sx ■ 


Proof For all x in cl(5a;), ^(x) is right invertible and we have ^(r*(x))(x) = 0. 


Thus, the rows of ^(t*(x)) are orthogonal to the column vectors of ^^(x). Moreover, they are 
independent due to the fact that it is a submersion. The Jacobian of t* can therefore be completed 
with §|-^(T*(a:)). □ 

Remark 2. Since is of constant rank n on O, the existence of such a function F is 
guaranteed at least locally by the constant rank Theorem. 

Example 2 (Continuation of Example]^. Elimination of the Xi in the 4 equations given by 
the injective immersion t* defined in ( |1.3[ ) leads to the function F(^) = ^ 2^3 satisfying (2.4). 

It follows that a candidate for completing: 


dr* 

dx 


(x) = 


( 1 
0 


V 


-xz 

0 


0 

1 

0 

-xz 


\ 


0 
0 

-x\ 
-X2 ) 


(2.5) 


IS 


dF 


(r*(x))^ = {x 2 Xz,-XiXz,X 2 ,-Xi)'^ 


Notice that this vector is nothing but the column of the minors of the Jacobian (2.5). It gives as 
determinant 


(0:2X3)^ + (xixz)'^ + xl+ xl 


which is never zero on O. 


With Proposition 2.1 
the function 


we know that, for any bounded open set Sx such that A C cl(iSa;) C O 


T*(a:, w) = (xi + X 2 Xzw, X 2 — XiXzw, —XiXz + X 2 W, —X 2 Xz — Xiw) 









is a diffeomorphism on Sx x Be{0) for e in ]0,e<>]. 

With this T*, (1.12) gives us the following observer written in the original a;-coordinates com¬ 
plemented with w : 


/ Xi 



1 

X3W 

X2W 

X2X3 

\ 

X3 



-X3W 

1 

—XiW 

-X1X3 


X2 



-X3 

w 

-Xl 

X2 


y w 



\ —w 

-X3 

-X2 

-Xl 

/ 


/ X 2 - X1X3W \ 


( ^ \ 


( ^ki\ 


—X1X3 + X2W 

—X2X3 — XiW 

V 0 ^ 

-f 

0 

0 

ysat(ii:E|)y 

-f 

ek2 

£^k 3 

——- 1 

1 


Unfortunately the matrix to be inverted is non singular for {x,w) in Sx x Be{0) only and we have 
no guarantee that the solutions of this observer remain in this set. This shows that another modi¬ 
fication of T* is needed to make sure that belongs to this set whatever ^ in is, namely 

to satisfy P2. A 


The drawback of this Jacobian completion method is that it asks for the knowledge of the 
function F. It would be nicer to have simply a universal formula relating the entries of the columns 
to be added to those of ^r. 


2.3. The P[m,n] problem. Finding a universal formula for the Jacobian completion comes 
back to solving the following problem, which we will call P[m, n] : 

For a continuous m x n matrix tp = {pij) of rank n (with n < m), can we find m — n continuous 
vector functions of the coefficients of ip, such that, with 'y{p) = ( 71 , ■.., 7m-n)) the 

matrix : 




( 2 . 6 ) 


is invertible ? 

It turns out that the existence of a solution to the P[m, n] problem is related to the simpler P[m, n] 
problem studied by Eckmann (among others). This problem is: 

For a continuous m x n matrix p = {pij) of rank n (with n < m), can we find a continuous vector 
function 71 of the coefficients of p which is orthogonal to each column of p ? 

Theorem 2.3. Eckmann theorem] The P[m, n\ problem is solvable if and only if the pair 
(to, n) is in the following table. 


m = 

even 

> 2 

7 

8 

n = 

I 

TO — 1 

2 

3 


The two most ’’common” cases are : 

• m = n + 1 : a solution is to add the vector of the corresponding minors (see Example 

• n = 1, TO even : a solution is to take permutations of the elements from the given vector. 
According to Theorem 2.3 the pairs for which the P\rn, n\ problem admits solutions are rather 

rare. This means that the existence of solutions to our P[to, n] problem is even rarer. Indeed, the 
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P[m, n] problem consists in adding only one column to a full-rank matrix, while we need m — n to 
solve the P[m, n\ problem. Thus, only the case m = n + \ solves our problem yet. 

In fact, when n = 1, i-e we need to find m — 1 vectors to complete the given one into an 
invertible matrix, looking for a universal formula comes back to answering the following question : 
For m > 3, do there exist m — 1 continuous vector functions 71 ,..., 7^-1 defined on K™ such that 
the matrix (z, 7 i(z),..., 7 m-i(-z)) is invertible for all z in \ {0} ? 

Restricting z to the unit sphere and projecting all the vectors "fiiz) onto the orthogonal complement 
of z, this is equivalent to asking whether the sphere S'^~^ is parallelizable (since the will be a 
basis for the tangent space at z for each z G It turns out that this problem admits solutions 

only in dimension m = 4 or to = 8 (see |12jb For instance, for to = 4, the matrix 


( 

Zi 

-^2 

2:3 

Zi \ 


Z2 

Zl 

-Zi 

zn 


z?, 

-Zi 

-Zl 

-Z2 

V 

Z4 

23 

Z2 

-Zl / 


is orthogonal for any (zi, Z2, z:3, Z4) in \ { 0 }. 

Noticing that the existence of solutions for P[to, n] implies the existence of solutions for P[to — 
1, n — 1], we deduce that neither P[7, 2] nor P[8,3] admits any solution. We finally conclude that : 

Theorem 2.4. The P[to, n] problem is solvable if and only if the pair (to, n) is in the following 
table. 


m = 

4 

8 

> 2 

n = 

1 

1 

TO — 1 


Fortunately, as illustrated in the following example, we can sometimes ’’modify” to and n in 
a way to get one of these good pairs even when it is not directly the case. Besides, we will see 
in Section how to change the initial immersion r* to have a universal explicit solution to the 
Jacobian completion problem. 

Example 3 (Continuation of Example . In Example we have completed the Jacobian 
(2.5) with the gradient of a submersion and observed that the components of this gradient are 
actually minors. We now know that this is consistent with the case to = n -I- 1. But we can also 
take advantage from the upper triangularity of the Jacobian (2.5) and complete only the vector 
{—Xi,—X 2 ) by for instance (x 2 ,—xi). Actually in doing so, we move to the case with n = 1 and 
TO even. The corresponding vector 7 is 7(0;) = (0, 0, X 2 , —Xi). Here again, with Proposition 2.1 we 
know that, for any bounded open set such that A C cl(5a;) C O the function the function 


Tf{x,w) = (Xi , X2 , -X1X3 + X2W , -X2X3 - Xiw) 

is a diffeomorphism on x Pe(0), where in this particular case s can be arbitrary, no need for it 
to be small. However, the singularity at x'l = X 2 = 0 remains and P2 is still not satisfied. A 

2.4. Wazewski theorem. Historically, the problem of Jacobian completion was first ad¬ 
dressed by Wazewski (see [22]). His formulation was : 

Given mn continuous functions (pij : O C K” —)■ K, look for m(m — n) continuous functions 
')ki : CJ —>■ M such that the following matrix is invertible for all x in O : 

P(x) = (vj(x) 7 (x)) . 


10 


(2.9) 











The difference with the previous section, is that here, we look for continuous functions of x instead 
of continuous functions of the components of ^p{x). This other version of the problem admits a far 
more general solution : 

Theorem 2.5 (|2H Theorems 1 and 3] and [71 page 127]). If O, equipped with the subspace 
topology of K”, is a contractible space, then there exists a C°° function 7 making invertible the 


matrix P{x) in (2.9) for all x in O. 


Remark 3. A key assumption in Theorem 2.5 is that the set O, on which the immersion t* is 
defined, is contractible. This is a strong requirement which for instance does not apply to any open 
set O contained in \ ({( 0 , 0 )} x M) containing the closure of A defined in (1.8) of Example^ 
We will see in Section^how we can change the immersion t* given by Assumptions A.II to satisfy 
this constraint. 

Proof. The reader is referred to 13 page 127] or [3 pages 406-407] and to (23 Theorems 1 and 
3] for the complete proof of existence of a continuous function 7 . Here are some sketches. 

The proof given by Eckmann in [7] can be summarized as follows. Let Vm,n be the space 
of orthogonal matrices of m rows and n columns, and let p be the projection which removes the 
m — n last columns of an orthogonal matrix of dimension m. It is known that {SO{m),p,Vm,n) 
is a Hurewicz fibration and thus we can ’Tift” the application ip : X ^ Vm,n to an application 
P : X —)■ SO{m) since X is contractible. Unfortunately, we have not been able to extract a 
“formula” from this proof. 

The proof given by Wazewski in |22] for spaces homeomorphic to a parallelepiped is more 
“workable”. It is based on the fact that, if ip{x) is a full rank m x n matrix, then the completed 
matrix (ip(x) jix)) is invertible whenever j(x) is a full-rank m x (m — n) matrix satisfying 


ipfxY"^{x) = 0 . 

To exploit this remark, we note that, if we have the decomposition 

A{x) 


( 2 . 10 ) 


ip{x) = 


B{x) 


with A{x) invertible on some given subset IR of O, then 


7 (x) = 


C{x) 

D{x) 


satisfies the above properties on 3? if and only if D{x) is invertible on 3? and we have 

C{x) = —{A^{x))~^B{x)^D{x) Vx € ?fi . 

Thus, C is imposed by the choice of D and choosing D invertible is enough to build 7 on 3?. 
Also, if we already have a candidate 


P{x) = 


A{x) Co{x) 
B{x) Dq\x) 


on a boundary 93R of 3?, then, necessarily, if A{x) is invertible for all x in 93R, then Dq{x) is invertible 
and Co (a;) = —(AA'{x))~^B{x)'^Dq{x) all x in diR. Thus, to extend the construction of a continuous 
function 7 inside 3R from its knowledge on the boundary 93R, it suffices to ensure D = Dq on OR. 
Because we can propagate continuously 7 from one boundary to the other, Wazewski deduces from 
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these two observations that, it is sufficient to partition the set O into adjacent sets where a given 
nxn minor Ai is invertible. This is possible since (f is full-rank on O. When O is a parallelepiped, he 
shows that there exists an ordering of the such that the continuity of each Di can be successively 
ensured. We illustrate this proof in Example below. 

It remains to complete the proof of Eckmann or Wazewski by showing how the continuous 
function 7 making P invertible can be modified into a smoother one giving the same invertibility 
property. Let 7 ^ denote the ith cofumn of 7 . 

We start with modifying 71 into 71 . Since ip, 7 and the determinant are continuous, for any x 
in O, there exists a strictly positive real number r^, such that 


det {(p{y) 71 (x) 72 :m-n(?/)) > 0 , Vy G Br,,(x) , (2.11) 

where 7 ^,^ denotes the matrix composed of the to columns of 7 . The family of sets 
is an open cover of O. Therefore, by |T31 Theorem 2.1], there exists a subordinate 
partition of unity, i.e. there exist a family of C°° functions ipx '■ O ^ K>o such that 

Supp C Sr, (a;) Vx G C> , (2.12) 

{Supp {iIjx)}x&o is locally finite , (2-13) 

My) = 1 Vy G O . (2.14) 

xGO 

With this, we define the function 71 on O by 

7i(y) = 51 . 

xeo 


This function is well-defined and C°° on O because the sum is finite at each point according to 


(2.13). Using multi-linearity of the determinant, we have, for all y in O, 


det {(f{y) 7 i(y) 72:m-„(y)) = 5 Z tp^(y) det {if{y) 71 (x) 72:m-„(j/)) 


X^Sx 


Thanks to (2.13), at each point y in Sx, there is a finite number of 'ipxiy) which are not zero. 


Moreover, according to (2.14), there is at least one ipxiy) strictly positive. On the other hand, with 
( 2 . 12 ) and ( 2 . 11 ), we know the determinant corresponding to any non zero ipx{y) is strictly positive. 
Therefore, we can replace the continuous function 71 by the C°° function 71 as a first column of 7 . 
Then we follow exactly the same procedure for 72 with this modified 7 . By proceeding this way, 
one column after the other, we get our result. □ 

Example 4. Consider the function 


/ 


ip{x) = 


1 

0 

- 2:3 

0 


VISN-51 gN-tl 


9p 


0 

1 

0 

-2:3 


0 

0 

-Xi 

-X2 

P 


\ 


/ 


where is a real number and p(xi,X 2 ) = max{0,1 — (x^ -I- x^)}^. p(x) has full rank 3 for any 
X in K^, since p(xi,X2) ^ 0 when xi = X2 = 0. Hence a continuous function exists to augment tp 
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into an invertible matrix. To follow Wazewski construction, let <5 be a strictly positive real number 
and consider the following 5 regions of 


= ] — oo, —d] X 
sftg = [-<5,5]2 X K , 
3?5 


3?2 = <^] X [<^) + 00 ] X K, 

= [—(5,<5] X [— 00 ,—d] X K, 

[( 5 , +oo[x]R^. 


We select S sufficiently small in such a way that p is not 0 in 3 ? 3 . 

We start Wazewski’s algorithm in lft 3 . There the matrix A is given by rows 1, 2 and 5 of and 
B by rows 3 and 4. With picking D as the identity, C is D gives rows 3 and 4 of 7 and 

C gives rows 1, 2 and 5. 

Then we move to the region 3 ? 2 - There the matrix A is given by rows 1, 2 and A oi p, B by rows 
3 and 5. Also £>, along the boundary between 3?3 and 3 ? 2 , is given by rows 3 and 5 of 7 obtained 
in the previous step. We extrapolate this inside 3^2 by kipping D constant in planes Xi =constant. 
An expression for C and therefore for 7 follows. 

We do exactly the same thing for 3 ^ 4 . 

Then we move to the region 3fii. There the matrix A is given by rows 1, 2 and 3 of <p, B by rows 4 
and 5. Also L>, along the boundary between 3fJi and 3 ^ 2 , between 3fii and 3^3 and between 3fii and 
3 ^ 4 , is given by rows 4 and 5 of 7 obtained in the previous steps. We extrapolate this inside 3fii by 
kipping D constant in planes X 2 =constant. An expression for C and therefore for 7 follows. 

We do exactly the same thing for 3 ^ 5 . 

Note that this construction produces a continuous 7 , but we could have extrapolated I? in a 
smoother way to obtain 7 as smooth as necessary. A 


3. About P2 : Image extension of a diffeomorphism. 


3.1. Problem definition and result. To satisfy condition P2 of Proposition o we may 
need to extend a given diffeomorphism r*, satisfying condition PI with some set Oe^ in such a way 
that the image of the extended diffeomorphism covers the entire K™. In order to keep PI, this new 
diffeomorphism must still be an extension of r* in the sense of ( 1 . 111 . 

There is a rich literature reporting very advanced results on diffeomorphism extension problem. 
In the following some of the techniques are inspired from [131 Chapter 8 ] and [131 pages 2, 7 to 14 
and 16 to 18](among others). Note however that we are interested in a weak version of this problem 
since we allow a small deformation of the set on which the given diffeomorphism is to be matched. 
In the following, we give a solution to this problem when the following property is satisfied. 

Definition 3.1 (Conditions B). An open subset E o/K"* is said to verify condition^ if there 
exist a function k : —>■ K, a vector field y, and a closed set Kq contained in E 

such that: 


1. E = {z& K"',k(z) < 0 } 

2 . Kq is globally attractive for y 

3. we have the following transversality assumption: 


F)k 


Vz G K™ : k{z) = 0 . 


^If not replace x by 


X 

\/i+RF‘ 
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Theorem 3.2 (Image extension). Let ij)' 'D C —>■ i’i'D) C K™ 6e a diffeomorphism. If 

ififD) verifies condition M or V is C'^-diffeomorphic to and if is C'^, then for any compact set 
K in V there exists a diffeomorphism ife : V ^ K"* satisfying : 


M-D) = 


lfe{z) = lf{z) Vz S K. 


From this theorem, we deduce the following corollary : 

Corollary 3.3. Let A he a bounded subset o/K”, Oe an open subset o/K™ containing Ax {0} 
and r* : Oe — t Tf{Oe) be a diffeomorphism satisfying (1.11) and such that either Tf{Oe) verifies 
condition B or Oe is -diffeomorphic to K"* and r* is C^. Then, there exists rf : Oe — 
that Tf satisfies Properties PI and P2. 

Remark 4. For Oe to be diffeomorphic to Oe must be contractible. 


\ such 


3.2. Proof of Theorem |3.2[ The proof of Theorem 3.2 uses the following two technical 
lemmas : 

Lemma 3.4 (simplified). Let E be an open strict subset of verifying B. Then, for any 
compact set K in E, there exists a diffeomorphism fi: K™ —)■ E, such that (j) is the identity function 
on K. 

See Appendixfor a proof with a more precise statement. In particular, the set K where is 
to be the identity function may not be compact. Besides, can be made as smooth as the vector 
field X iR condition B. This proof is partly constructive as illustrated in the following example. See 
also Section [5l 

Example 5. Let M be an m x m matrix. Consider the function F : K 

and let the sets E and K be 

E = <5] , a:= M™,E(0^ = 0} 

where 5 is some strictly positive real number. E verifies condition B with 

T 




x(6 = -i ■ 


Indeed, the origin is in E and is asymptotically stable for y. Moreover, on the boundary of E, 
where 

the transversality assumption is verified since we have 

-^iOxiO = -s<o. 

To get a diffeomorphism, first, we define a ’’layer” close to the boundary of E as 




1)4 -4e . 


^ ) < (5 


It is the set of all the points reached at some time t in (0, e) by a solution of 

? = x(C) 
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(3.1) 




issued from a point in the boundary of E. It is contained in E. We want the diffeomorphism to be 
the identity function on K. For this we consider the complementary set within E of the “layer” 


= U e 




It contains K. To construct the diffeomorphism, we associate to any point ^ in \ E^ 


4 5 


tx is the time a solution of (3.1) issued from ^ needs to reach the boundary of E. With this we 
define the diffeomorphism (j) : R™ E a,s : 

m = 

where ^ : R - 




if < e-4"<5, 


g *£ ^ otherwise, 

is any C^, strictly decreasing function satisfying : 
v(t) = —t yt < —e , 

For example it can be 


(3.2) 


lim v(t) = 0 . 

t^OO 


v{t) = 


2 £ -|- t 


Wt>—e. 


(3.3) 

A 


Xi 



Fig. 3 . 1 . Construction of the diffeomorphism in Example^ 

Lemma 3.5. Consider a diffeomorphism ip : Br{0) — t ip{Bii{0)) C R™. For any e strictly 
positive, there exists a diffeomorphism iff, : R™ —>■ R™ satisfying 

1pe(z) = 1p{z) Vz G cI{Br-s{ 0 )) . 
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A shown in Appendix]^ this result is a direct consequence of [H Theorem 8.1.4], Unfortunately 
its proof is less constructive than for the previous case. 


Proof, of Theorem 3.2 


First case: ip(F) satisfies B; ij) being a diffeomorphism on an open set V, the image of any compact 


subset AT of I? is a compact subset of According to Lemma 3.4 there exists a diffeomorphism 

(p from to 'ipi'D) which is the identity on ip{K). Thus, the extension tjjg = cj)~^ o tp satisfies the 
required properties. 

Second case: T> is -diffeomorphic to K™ and ip is : Let (f>i : V ^ K"* denote the corresponding 
diffeomorphism. Let Ri be a strictly positive real number such that the open ball (0) contains 
4>i{K). Let i ?2 be a real number strictly larger than i?i. With Lemma 3.4 again, and since 
verifies condition B, there exists of -diffeomorphism (j )2 '■ Br^{Q) —>■ K™ satisfying 

(j) 2 {z) = z Vze Br^{0) . 


At this point, we have obtained a C^-diffeomorphism (p = <p 2 ^ ° <pi : T) ^ Br^{0). Consider 


X = P;o<p-^ :BrM _ 

that Xe = tp o on Br^ (0). Finally, consider ipe = Xe o (pi : V 


ippV). According to Lemma 3.5 we can extend A to Ag : M’” —>■ K™ such 

R™. Since, by construction of 


02, <^ = 01 on (pi (S_Ri(0)) which contains K, we have 0e = "0 on K. □ 


3.3. Relaxing con dition B. The representation of E in terms of the function k, in con¬ 
dition B for Lemma 3.4 can be replaced by the following weaker condition B’'®0 

Definition 3.6 (Relaxed condition B’’'^*). An open bounded subset E o/K"* is said to verify 
the relaxed condition if there exist a bounded vector field x, o,n-d a compact set Kq contained 
in E such that: 

1. Kq is globally asymptotically stable for x 

2. E is forward invariant by x- 

The following lemma shows that the relaxed condition B""®* implies the condition B up to a 
small deformation of E : 

Lemma 3.7. Let E be an open bounded subset o/K™ verifying the relaxed condition For 
any strictly positive real number d, there exists a bounded set E such that 


:1{E) C £ C {z e M™, d{z,E) G [0,d]} 


and E verifies condition B. 

Proof. The compact Kq being globally asymptotically attractive and interior to E which is 
forward invariant, E is globally attractive. It is also stable due to the continuity of solutions with 
respect to initial conditions uniformly on compact time subsets of the domain of definition. So it 
is globally asymptotically stable. It follows from [23l Theorem 3.2] that there exist C°° functions 
Vk : R™ —>■ R>o and Ve ■ R™ —>■ R>o which are proper on R™ and a class /Coo function a satisfying 

a{d{z,Ko)) <Vk{z) ,VKiz) = 0 W z G Kq , 
a{d{z,E)) <Ve{z) ,Ve{z) =0 V z G E , 

^{z)x{z)<-Vk{z) VzeR"^ 

^{z)x{z)<-Ve{z) VzGR™. 
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With the notations ve = E)<d^K{z) and fi = and since ot is of class K,oo^ we 

obtain the implications 

VEiz)+fiVKiz) = a(d) a{d{z,E))<VE{z)<a(d) , 

^ d{z, E) < d , 

^ VKiz) < VE ■ 


With our definition of /i, this yields also 

a(d) -^iVk{z) = Ve{z) 0 < ^< Ve{z) , 

^ 0 < d{z, E) < d . 


On the other hand, with the compact notation V{z) = Ve{z)+^Vk{z), we have ff (z) x(z) < —V{z), 
for all z S M™. All this implies that the sublevel set £ = {z € : V{z) < a{d)} is contained in 

{z G K’" : d{z,E) G [0, d]} and that cl(_E) is contained in £. Besides, £ verifies condition B with 
the vector field y and the function k = V — a(d). □ 

Suppose we have a diffeomorphism t* : ^ T*{Oe) where T*{Oe) is a strict subset of R™. To 

satisfy property P2, we need to modify this diffeomorphism. With Corollary |3.3[ Lemma [3.4| and 
this is possible if T*{Oe) verifies condition B or Og is C^-diffeomorphic to 


Lemma 


3.7 


and 


is C^. But for this, the given diffeomorphism r* needs to be defined on an outer approximation of 
Oe- If this is not the case, we can take advantage of the following lemma : 

Lemma 3.8. Let a diffeomorphism if be defined on an open set E m R^. Consider a bounded 
open set E' such that cl[E') C E. There exists a strictly positive real number dg such that for all 
X in E', 


d[if{x),W^\if{E))>do 


Proof. Suppose that, for each integer k, there exists Xk in E' such that 

d {il}{xk),M.'^\if{E)) < i . 

By compactness, there exists x* in cl(i?') C E such that d {ip{x*),'R"^\ijj{E)) = 0. Thus, ip{x*) is 
in cl{M.'^\ijj{E)) C]if{E). But since E is open, if{E) is also open by Brouwer’s invariance theorem. 
Thus, cl(R"*\'(/;(A)) n ip{E) is the empty set and we have a contradiction. □ 

To round the difficulty mentioned before this Lemma, we consider a bounded open set O' 
satisfying A x {0} C cl(O') C Oe- According to Lemmathere exists dg such that, 

d{ix,w),R^\{Oe))>do , V(cr,w;)eO', 

d{r:{x,zc),(r: (R™\(Oe))) > dg , y{x,w) G O' . 


Thus, i f Tf {Oe) verifies condition B or Oe is O^-diffeomorphic to R"* and t* is O^, by applying 
Lemma 3.7 with d < dg to r* (Oe) or to Oe, we get the existence of an open bounded set Of verifying 
A X {0} C cl(Og) C Og C Oe and such that TfjO l) verifies condition B or Og is O^-diffeomorphic 
to R™. We are now able to apply Corollary 3.3 and obtain a diffeomorphism r* : Oi —)■ R™ the 


We are now able to apply Corollary 3.3 
image of which is R"*. 
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3.4. DifFeomorphism extension from a snbmersion. As in Section 2.4 we consider now 
the particular case where there exits a full-rank function F : K™ —>■ K"*”" such that 

F{t*{x)) = 0 Wx € Sx ■ 


Thanks to Propositions 2.1 and |2.2[ for e sufficiently small, the function 

f)P 

(x,w) !->■ T*(x,w) = T*(x) -\- —{t*(x))'^W 


(3.4) 


is a diffeomorphism on x B^{0) but the set t*{Sx x Se(0)) may not be the entire K’" as required 
for P2 to hold. 

Proposition 3.9. Assume the existence of a full-rank function F : K™ — > K'"”" such that 
the set 




(3.5) 


is compact and we have 


T*{Sx) = {f€W^, F(C)=0} . 


(3.6) 


Let also e be such that the function defined in j[3.4\ ) is a diffeomorphism on Sx x Under these 

conditions, there exists a strictly positive real number d such that the set 

= < < 5 } 


satisfies 


T*iSx)CE^ CrfiSxXBM)- (3.7) 

Moreover if F is on K™, then E^ verifies condition B. 

Proof. We start by showing the existence of a strictly positive real number r such that 

e K™ : d{f,T*{Sx)) <r}c TfiSx x BM) ■ (3-8) 

We know that for all in T*iBx), there exists Xq in Sx satisfying fo = r*(a;o,0). Moreover, the 
Jacobian of r* at (xq, 0) is invertible. It follows from the Implicit function Theorem that, for each 
fo in T*{Sx), there exists a strictly positive real number rj such that, for all f in t*{Sx), satisfying 
1^ — Col < there exists {x,w) in Sx x Be{0) satisfying Tf{x,w) = f. Existence of r follows from 
compactness of cl(iS 3 ,) and therefore of cl(T*(iSa;)). 

Secondly, there exists J > 0 such that E^ is contained in {C G : d{f,T*{Sx)) < r}. Indeed, 
assume this is not the case. Then, for each integer k, there exists Cfc satisfying 


\Eifk)\<^ , d{^k,T*{Sx))>r. 


The sequence being in the compact set defined in (3.5), there is a subsequence which converges to 
C* satisfying 


d{e,r*{Sx))>r. 


F{C) = 0 
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This contradicts (3.6). Then (3.7) follows from (3.6) and (3.8). 


Finally, for the last point, we note that, with x the vector field defined by 


f9F OF ^ 


-1 


m, 


we have, for all ^ in K™, (^F^F') (^) = —F{^)^F{^). This implies that the closed set t*{Sx) is 

globally attractive for x and that the transversality assumption is verified. Thus, satisfies B. □ 
The example below illustrates the fact that, if we have a submersion F giving an (outer) 
approximation only of t* (Sx x Sw) we can still apply the diffeomorphism image extension of Theorem 

[Q 

Example 6 (Continuation of Example]^. In Example]^ we have introduced the function 

F{0 = 66 - 66 = 

as a submersion on M^\{0} satisfying 


F(r*(6)=0. 


(3.9) 


With it we have “augmented” the injective immersion t* given in (1.3) as 


dF^ 

T*{x, w) = T*{x) + {t*{x)) W = T*{x) + Mt*{x) W 

which is a diffeomorphism on Oe = SxX] — s^, —£<>[ for some strictly positive real number — £«. To 
satisfy P2 we need an image extension. 

It is difficult to obtain an exact representation of the image T*{Oe), but we can use F to 
obtain an approximation. Indeed, with (3.9), we obtain F[t*{x,w)) = \t*{x)\'^w or equivalently 
w = ^on zero, i.e. away from the observability singularity xi = X 2 = 0. 

Hence ensuring that ^ = t*{x,w) remains in 

= {^G K’”,F(6^ < <5} 

ensures that w remains small unless x is close to the observability singularity. 

In Example]^ we have obtained in (3.2) an expression of a diffeomorphism —>■ which 

is the identity on 


With it, we can replace t* : Oe —)■ T*{Oe) by (j)~^ or*, ensuring that, for any {x^w) in t*{x,w) 
is in E^ for S sufficiently small. In doing so, the domain of non invertibility of the Jacobian has 
been reduced. But we still have a problem for (xi,X 2 ) close to the origin where the observability 
singularity occurs. A 
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4. Modifying A.II and A.C to allow PI and P2 to be satisfied. Throughout Sections 
[^and|^ we have studied tools which enable us to fulfill conditions PI and P2 of Proposition |1.1| 
However, they need conditions on the dimensions or on the domain of injectivity O which are 
not always satisfied : contractibility for Jacobian extension and diffeomorphism extension, limited 
number of pairs (m, n) for the P[to, n] problem, etc. Expressed in terms of our initial problem, 
these conditions are limitations on the data /, h and r* we considered. In the following, we show by 
means of examples that, in some cases, these data can me modified in such a way that our various 
tools apply and give a satisfactory solution. Such modifications are possible since we restrict our 
attention to system solutions which remain in A. Therefore we can modify arbitrarily the data /, h 
and T* outside this set. For example we can add arbitrary “fictitious” components to the measured 
output y as long as they are constant on A. 

4.1. For contractibility. It may happen that the set O attached to t* is not contractible, for 
example due to an observability singularity. We have seen Jacobian completion and diffeomorphism 
image extension may be prevented by this. A possible approach to round this problem when we 
know the system solutions stay away from singularities is to add a fictitious output and to modify 
Assumptions A.II and A.C. 

Example 7 (Continuation of Example]^. The observer we have obtained at the end of 
Example for the harmonic oscillator with unknown frequency is not satisfactory in particular 
because of the singularity at xi = X 2 = 0. To round this difficulty we add, to the given measurement 
y = xi, the following “fictitious” one : 


y2 = h2{x) = p{xi,X2)[x:i-<,] 

with in ]0,r[ and p{xi,X 2 ) = maxjO, ^ — {x^ +x|)}^. By construction this measurement is 
zero on A meaning that we do not change the data on this set. The interest of y 2 is to give us 
access to X 3 directly. Indeed, consider the new function t* defined as 


T*{x) = (xi , X2 , -X1X3 , -X2X3 , p(xi,X2) [X3 - <r]) . 
T* is on and its Jacobian is : 


/I 0 0 \ 

0 1 0 

—X 3 0 —xi 

0 -X3 -X2 

P y 


(4.1) 


(4.2) 


which has full rank 3 on K^, since p(xi,X 2 ) ^ 0 when xi = X 2 = 0. It follows that the singularity 
has disappeared and this new r* is an injective immersion on the entire which is contractible. 

We have shown in Example how Wazewski’s algorithm allows us to get in this case a 
function 7 : —>■ satisfying : 


det 




= 0 


Vx G . 


This gives us t*{x,w) = t*(x) + 'y{x)w which is a C^-diffeomorphism on x with e 

sufficiently small. _ 

Furthermore, Og = x Be{0) being diffeomorphic to K®, Corollary 3.3 applies and provides a 
modification r* of r* satisfying PI and P2. A 
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4.2. For a solvable P[m,n] problem. As already seen in Example]^ it is appropriate to 
exploit the fact that some rows of , say p of them, are independent for all x in A. Indeed this 
allows to transform the Jacobian completion problem into a P[m — p^n — p] problem. If the new 
pair {m — p,n — p) is not in table (|2.8[) but {m — p + k,n — p), we can still add k arbitrary rows to 


dr 

dx 


Example 8 (Continuation of Example]^. In Example]^ by adding the fictitious measured 
output 2/2 = ^2 (x), we have obtained another function r* for the harmonic oscillator with unknown 
frequency which is an injective immersion on In this case, we have n = 3 and m = 5 which 
gives a pair not i n ta ble (2.8). But, as already exploited in Example]^ the first 2 rows of the 
Jacobian in (4.2) are independent for all x in It follows that our Jacobian completion 
problem reduces to complement the vector 


{-Xi,-X2,p{Xi,X2)) 


This is a ^*[3,1] problem which is not in the table (2.8). Instead, the pair (4,1) is in the table, 
meaning that the following vector can be completed via a universal formula 


(-a;i,-a;2,p(a:i,X2),0) . 

We have added a zero component, without changing the full rank property. Actually this vector is 
extracted from the Jacobian of 


T*{x) = {xi , X 2 , -xiX:i , -X 2 X 3 , , p(xi, X 2 ) [ 2:3 - <?], 0) (4.3) 

where the added last component at 0 is consistent with the high gain observer paradigm, when we 
add another fictitious measured output 


ya = 0 . 

A complement is 

/ 2:2 -p 0 \ 

-2;i 0 -p 

0 —a;i —X2 

\ p 2:2 -2;i / 


It gives the diffeomorphism 

tI{x,w) = (^Xi, X2, [-2:ia;3 + X 2 W 1 - p{xi,X2)W2] , [- 2 : 2 X 3 - XiWi - p(Xi,X2)w3] , 

[p(xi,X 2 )[x 3 - - X 1 W 2 - X2'U:3] , [p(xi,X 2 )'u:i + X 2 W 2 - XiWs)]^ . 

the Jacobian of which is : 


drl 

dx 


(x,w) = 


1 0 

0 I 

dlD 

- ^W2 

-^1 - Ia^3 -X3 - 

^[X3 -i]-W2 ||[X3 -<i]-W3 

V - ^3 + W 2 
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0 0 0 0 \ 
0 0 0 0 
-Xi X 2 -p 0 

-X2 -xi 0 -p 

P 0 — Xi — X2 
0 P X2 -XiJ 









with determinant 


(x^+xl + p(xi,X2ff 

which is nowhere 0 on K®. So r* is a diffeomorphism defined on K® which is also surjective since it 
is proper (Hadamard-Levy Theorem) thanks to 

xi =ii , X 2 = i 2 , 


/ 

6 

-p{iiii2) 

0 

\ 

/ 2:3 \ 


/ 6 \ 

-?2 


0 

-p(li,6) 


Wi 


I 4 

P(ll,l 2 ) 

0 


-I 2 


W2 


Is 

V 0 

P(|li| 2 ) 

6 

-6 

/ 

\W 3 J 


\ le / 


So, with the addition of the fictitious measured outputs 2/2 = ^ 2 ( 2 ^) and 2/3 = 0, we have obtained 
the new function r*, given in ( |4.3[ ), to be used, in place of (1.31, as another starting point for the 
construction of an observer and in particular for the construction of a function t* satisfying the 
properties PI and P2 of Proposition [T]^ Also to obtain a set satisfying the property P3, we can 
follow the high gain observer paradigm and complete the expression of ip in taking advantage 
of the fact that, for x in A, we have 


h = p{xi,X2)[x^ - c] = 0 

This motivates the new function p 

( 


2/3 = 0 


Vitx,y) = 


i2+£ki{y-xi) \ 
is + ^‘^k2iy - xi) 
ii + i^k3{y - xi) 
sat(xii§) + l‘^ki{y — xi) 

—a is 

V -H. 


where the function sat is defined in (1.9) and a and b are arbitrary strictly positive real numbers. 
With picking I large enough, it can be paired with any function r : K® —K® which is locally 
Lipschitz. With all this, we have obtained the following observer for the harmonic oscillator with 
unknown frequency 


/ 


X 2 


0 

1 

0 

0 

0 

0 

£3 


-*3 - §f^W 2 

- if 

-xi 

£2 

-p 

0 

Wi 


-^1 - if ^3 

-*3 - §^,W3 

-£2 

-£i 

0 

-p 

W2 


§^[X3-<^]-W2 

if[^3-?]-W3 

P 

0 

-£i 

-£2 

[^ 3 / 


CO 

1 

,^|cg 

if^l +^2 

0 

p 

£2 

-£ij 


X 


(4.4) 


/ 


^2-r^fcl(2/-il) \ 

[-X1X3 + * 2^1 - pixi,X2)w2] + £'^k2{y - Xi) 
[-^ 2*3 - XiWi - p{xi,X2)w3] + tk3{y - ii) 

sat(a;ii|) + i'^k4{y — ii) 

-a [p{xiA2)[x3 - ?] - ^ 1^2 - ^ 2 ^ 3 ] 

, -b[p{xiA2)wi+X2W2- X1W3)] ) 
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It is globally defined and globally convergent for any solution of the oscillator initialized in 

A = {a; e : x\+ xl & ]l,r[ , S]0,r[} . 

As a final remark, we observe that the manifold Cs = Ce = 0 is invariant. This implies the 
existence of an observer with order reduced to 4. Unfortunately it cannot be expressed with 
coordinates c = {xi,X 2 ,x^,w) since the following Jacobian is singular for = ^2 = 0 : 


dc 


/ 1 0 0 0 
0 10 0 


where * denotes an unimportant term. 

The observer (4.4) is just an illustration of what can be obtained by using in a very nominal 
way our tools. We do not claim any property for it. For example, by using another design, an 
observer of dimension 2, globally convergent on A, can be obtained. 

A 


5. Example of the bioreactor. As another illustration we consider the model of bioreactor 
presented in [S] and modeled by the dynamics : 

. Q^\X\X2 . Qj'^Q.\X\X2 

Xi = - UXi , X2 = - UX2 + UQa , y = Xi 

a2Xi + X2 a2Xi + X 2 

where the a^’s are strictly positive real numbers and the control u verifies : 0 < Umin < 'a(t) < 
Umax < fli- This system evolves in the set O = {x € : Xi > £i , X2 > —02X1} which is forward 

invariant. 

A high gain observer design leads us to consider the function r* : CJ —?> defined as : 


T* {XI,X 2 ) 



aiXiX2 
a2Xi + X2 


It is a diffeomorphism onto 


T*{0) = {eeM' : Cl >0, >6} ■ 

The image by r* of the bioreactor dynamics is of the form 

ii = 6 +5i(6)w 
6 = ‘^771(^1,6) + 52 (Ci)C 2 )m 


for which the following high gain observer can be built (see [5]): 

| = + (5.1) 

where C = (liO)"^ and S^o is solution of AA'Soo + S^oA — C^C = —fS'oo for I sufficiently large. 
Since t* is a diffeomorphism, the expression of this observer in the cc-coordinates is : 


X 


a\X\X2 

a2Xi+X2 


— UXi 


a3aiXiX2 

a2Xi-\-X2 


— UX2 +ua4 


aia2xl 0 \ - xi)) 

—aia2x\ {a2Xi+X2Y J aia2x\ 
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Unfortunately the right hand side has singularities at aii = 0 and X 2 = —aiXi and nothing guaran¬ 
tees that O is still invariant for the observer dynamics (5.2). To round this difficulty, we can modify 
T* to make its image equal to while keeping it unchanged inside a subset of O, arbitrarily close 
to O. 

To construct this modification, we view the image t*{ 0 ) as the intersection t*{ 0 ) = Eif] E 2 
with : 


^1 — {(?l!C2) G > £ 1 } 


E2 = {(?i,C 2 )gK^, aiCi > C2} 


This exhibits the fact that t*{ 0 ) does not satisfy the condition B since its boundary is not C^. 
We could smoothen this boundary to remove its “corner”. But we prefer to exploit its particular 
“shape” and proceed as follows : 

• first, we build a diffeomorphism —?> Ei which acts on without changing ^ 2 - 

• then, we build a diffeomorphism (/)2 : —>■ E 2 which acts on ,^2 without changing ^ 1 . 

• finally, we take t* = o o t* : O — )■ 

To build (j)i and ^ 2 , we follow Lemma |3.4| since Ei and E 2 satisfy condition B with : 


Ki(0 = £1 - Cl, ^^ 2(0 = 6 - OiCi , Xi(C) = 


-(Ci-1) 

0 


, X2(C) = 


0 

-(C 2 + 1 ) 


By following the same steps as in Example]^ with e an arbitrary small strictly positive real number 
and V defined in (3.3), we obtain : 


^€,1 - 


Inip^ 

l—£ 


Ee,i — {(Cii C 2 ) G Cl > 1 ~ (1 ~ £)e 


^i(C) = 


c 


_ 1 ) + 1 


ifCe^e.i 

otherwise 


t^,2 = 


ai?l-l-l ’ 


(5.3) 


E, 


£,2 


= {(Cl,C 2 )GK^ C 2 <e-naiCi + l)-l} 


02(C) = 


c 


-f 1) - 1 


ifCGi^e.2 

otherwise 


We remind the reader that, in the ^-coordinates, the observer dynamics are the same let it be 
built with T* or with r*. The difference is only in the way x is related to Moreover this “way” 
is such that it has no effect on the system solution we have 

r*(x)=7^(x) Va;€C>“-£”. 

As a consequence the difference between t* and r* acts only during the transient, making sure that 
X never reaches a singula rity of the Jacobian of r*. 

the results in the ^ coordinates (to allow us to see the effects of both 


5.1 


We present in Figure 

T* and T*) of a simulation with (similar to [5]) : 

Ui = U2 — ^3 — 1 5 U4 = 0.1 
u{t) = 0.08 for t < 10 , = 0.02 for 10 < t < 20 , 

a;(0) = (0.04,0.07), x(0) = (0.03,0.09), 
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= 0.08 for t > 20 
i = 5. 








T rajectory in the image 



Fig. 5 . 1 . Bioreactor and observers solutions in the ^-coordinates 


The solid black curves are the singularities locus. The red (= solid dark) curve represents 
the bioreactor solution. The magenta (= light grey dashdot) curve repre sent s the solution of the 
observer built with t*. It evolves freely in according to the dynamics (5.1), not worried by any 
constraints. The blue (= dark dashed) curve represents its image by (j)~^ which brings it back inside 
the constrained domain where r = can be used. This means these two curves represent the 
same object but viewed via different coordinates. 

The solution of the observer built with r* coincides with the magenta (= light grey dashdot) 
curve up to the point it reaches one solid black curve of a singularity locus. At that point it leaves 
T* (O) and consequently stops to exist in the ii-coordinates. 


As proposed in [Mia, instead of keeping the raw dynamics ( |5.1[ ) untouched as above, we can 
modify them to force ^ to remain in the set t*{0). Taking advantage of the fact that this set is 


convex, the modification proposed in [3] consists in adding to (5.1) the term 




(5.4) 


max{Ki(^) + S, 0}^ 
max{K2(|) + S, 0}^ 

large real number. The solution corresponding to this modified observer dynamics is shown in 


where 1 )(^) = 


with 6 an arbitrary small real number and 7 a sufficiently 


Figure 5.1 with the dotted black curve. As expected it stays away from the the singularities locus 
in a very efficient way. But for this method to apply we have the restriction that r* (O) should be 
convex, instead of satisfying the less restrictive condition B. Moreover, to guarantee that ^ is in 
T*{0), 7 has to be large enough and even larger when the measurement noise is larger. On the 
contrary, when the observer is built with r*, there is no need to tune properly any parameter to 
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obtain convergence, at least theoretically. Nevertheless there maybe some numerical problems when 
^ becomes too large or equivalently is too close to the boundary of t*{0). To round this 

problem we can select the ’’thickness” of the layer (parameter e in (5.3)) sufficiently large. Actually 
instead of “opposing” the two methods, we suggest to combine them. The modification (5.4) makes 
sure f does not go too far, and t* makes sure that x stays away from the singularities locus. 


6. Conclusion. We have presented a method to express the dynamics of an observer in the 
given system original coordinates enlarging its domain of validity and possibly avoiding the difficult 
left-inversion of an injective immersion. 

It relies on the assumption that we know an injective immersion and a converging observer for 
the immersed system. The idea is not to modify this observer dynamics but to map it back to the 
original coordinates in a different way. Our construction involves two tools : the extension of an 
injective immersion into a diffeomorphism through Jacobian completion and the extension of the 
image of the obtained diffeomorphism to enlarge the domain of validity. 

For Jacobian completion, we are relying on the results of Wazewski [22] and Eckmann [7]. It 
allows us to build a diffeomorphism by complementing the given coordinates with new ones and to 
write the given observer dynamics in these augmented coordinates. 

For diffeomorphism extension, we have proposed our own method widely inspired from m 
Chapter 8] and m pages 2, 7 to 14 and 16 to 18]. 

In our presentation we have assumed the system is time-invariant and autonomous. Adding 
time-variations is not a problem but dealing with exogenous inputs is more complex. This is in part 
due to the fact that, as far as we know, the theory of observers, in presence of such inputs, relying 
on immersion into a space of larger dimension, as high gain observers or nonlinear Luenberger 
observers, is not available yet. Progress on this topic has to be made before trying to extend our 
results. 


Appendix A. Construction of a diffeomorphism from an open set to K'". In this 
section we prove a result slightly more general than Lemma 3.4 For this we use the following 
notations: 


The complementary, closure and boundary of a set S are denoted 5“ 
The Hausdorff distance dn between two sets A and B is defined by : 


cl(S') and dS, respectively. 


dH{A,B) = max < sup inf |z^ —zb|,sup inf \za — zb\\ 

With Z{z,t) we denote the (unique) solution, at time t, to i = xi^) going trough x at time 0. 

Lemma A.l. Let E be an open striet subset of ME' verifying B, with a C® veetor field y. Then, 
for any strictly positive real number e, there exists a C'^-diffeomorphism (f: K™ -A E, such that, 
with 

E= U Z{dE,t) , 

tG[0,£] 

we have (j){z) = z for all z & E^ = E and dH{dE^,dE) < esupjx(z)|. 

Proof. We start by establishing some properties. 

- E is forward invariant by x- This is a direct consequence of points and of the condition B. 

- E is closed. Take a sequence (zfc) of points in E converging to z*. By definition, there exists a 
sequence (t^), such that : 

tfc G [0,e] and Z{zk,—tk)&dE V/c G N . 
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Since [0,e] is compact, one can extract a subsequence (^^(fc)) converging to t* in [0,£], and by 
continuity of the function {z,t) !->■ Z{z,—t), {Z{z^i^k),ta(k))) tends to Z{z*,—t*) which is in dE, 
since dE is closed. Finally, because t* is in [0,e], z* is in S by definition. 

- Ti is contained in cl{E). Since, E is forward invariant by y, and so is cl(£') (see |lll Theorem 
16.3]). This implies 


dE c Z.= IJ Z{dE,t) c c1(F;) = EUdE. 

te[o,e] 

At this point, it is useful to note that, because E is a closed subset of the open set E, we have 
E n A = E\9A. This implies : 

E\E^ = {E^Y f^E = {E^l^T)f^E = Y.f^E = E\dE, (A.l) 

and E = E,U (E\5A). 

With all these properties at hand, we define now two functions tz and 9z- The assumptions of 
global attractiveness of the closed set Kq contained in E open, of transversality of x to dE, and the 
property of forward-invariance of E, imply that, for all z in E'^, there exists a unique non negative 
real number tz satisfying: 


K{Z{z,tz))=0 Z{z,tz)€dE. 

The same arguments in reverse time allow us to see that, for all z in E, tz exists, is unique and in 
[—£,0]. This way, the function z —)■ is defined on {E^Y. Next, for all z in {E^Y, we define : 

9z = Z{z,tz). 

Thanks to the transversality assumption, the Implicit Function Theorem implies the functions 
z tz and z 9z are C"* on {E^^Y■ 

Remark 5. n having constant rank 1 in a neighborhood of dE, this set is a closed, regular 
submanifold of K™. The arguments above show that z i—)■ (0^, tj,) is a diffeomorphism between E'^ 
and dE x [0, +oo[. Since dE is a deformation retract of E’^ and the open unit ball is diffeomorphic 
to K™ [To], if E were bounded, E^ could be seen as a h-cobordism between dE and the unit sphere 
and tz as a Morse function with no critical point in E’^. See [19] for instance. 

Now we evaluate tz for z in dT,. Let z be arbitrary in 9E and therefore in E which is closed. 
Assume its corresponding tz is in ] — e, 0[. The Implicit Function Theorem shows that z tz 
and z 9z are defined and continuous on a neighborhood of z. Therefore, there exists a strictly 
positive real number r satisfying 

Vy G Br{z) , 3ty g] - £, 0[ : Z{y, ty) G dE . 

This implies that the neighborhood Br{z) of z is contained in E, in contradiction with the fact that 
z is on the boundary of E. 

This shows that, for all z in 9E, tz is either 0 or —e. We write this as 9E = dE U (9E)j, with 
the notation ((?E). = |z G S : tz = 

Now we want to prove dE,, C (i9E)^. To obtain this result, we start by showing : 

dE^ fidE = tb and dE^ c 5E . 
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(A.2) 


Suppose the existence of z in dE;. n dE. z being in dE, its corresponding is 0. By the Implicit 
Function Theorem, there exists a strictly positive real number r such that, 

Vy e Br{z) , e ]-§, |[ : Z{y,ty) e dE . 


But, by definition, any y, for which there exists ty in ] ^, 0], is in S. If instead ty is strictly positive, 
then necessarily y is in E'^, because E is forward-invariant by y and a solution starting in E cannot 
reach dE in positive finite time. We have obtained : C E U £1^ = {E^y. Br{z) being a 

neighborhood of z, this contradicts the fact that z is in the boundary of E^. 

At this point, we have proved that dE^ n dE = 0, and, because E^ is contained in E, this 
implies dE^ C E. With this, (A.2) will be established by proving that we have dE^ C dE. Let z 
be arbitrary in dE^ and therefore in E which is open. There exists a strictly positive real number 
r such that we have : 


Br{z) nE^ = Br{z) n (£; n E=) 7 ^ 0 , b^z) n = Br{z) n (£;= u e) 7 ^ 0 , b^z) c e . 

This implies Br{z) n E° 7 ^ 0 and Br{z) n E 7 ^ 0 and therefore that z is in dE. 

We have established dE^ n dE = 0, dE^ C 9E and dE = dE U (9E)-. This does imply : 

dE^ C {dE)^ = {z&E -. t, = -e} . (A.3) 


This allows us to extend by continuity the definition of tz to K™ by letting tz = —e for all z G Eg. 
Thanks to all these preparatory steps, we are finally ready to define a function cj) : K™ —>■ E as: 


(/)(z) 


Z (z, tz + v{tz)) , if Z e {EgY , 
z, \i z G Eg , 


(A.4) 


where v is an arbitrary C® and strictly decreasing function defined on K satisfying: 


v{t) = —t Vt < —e , lim u(t) = 0. 

t—>- + oo 

The image of (j) is contained in E since we have Eg G E and : 


tz + y{tz) > tz 'izG El , 
Z{z,tz) G dE , 

Z{z,t) G E y{z,t) G dE X Ryo . 


The continuity of the functions (z,t) S K™ x K i-G Z{z,t) G M. and z G E^ ^ tz G [—£,-1-00 [ 
implies that this function (f) is continuous at least on W^\dEg. Also, for any z in dEg, tz is defined 
and equal to —e (see ( A.3[ )). So, for any strictly positive real number rj, there exists a real number 
r such that : 


\ty + e\ < 
v{ty)+e < 

(t>{y) = 
yiy) = 


V 

y 


z{y,ty + v{ty)) 


Vy e Br{z) , 

Vy G Br{z) , 

Vy G Briz) n Eg , 
yyGBr{z)GEl . 


Since we have : 


(j){z) = Z{z,tz + y{tz)) = Z {z,-e + v{-e)) = z. 
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we conclude that 4> is also continuous at z. 

By differentiating, we obtain : 

- at any interior point z of 

- at any z in Eg (which is open) ^(z:) = I. Also, for any z in dEg^ we have : 

rl7 r)f r)7 r)f 

-^{Z,E + HQ) + x(Z(z, t. + HQ))^(z)(1 + ^'(Q) = ^(z, 0) + x(Z(z, 0))^(z)(l - 1) , 

= / . 

This implies that (j) is Q on K"*. 

We now show that (p is invertible. Let y be arbitrary \n E 7 E^ = E 7E. There exists ty in 
[—£, 0[. The function v being strictly monotonic, v~^(ty) exists and is in [—e,+oo[. This allows us 
to define properly (p~^ as : 


/ Z{y,ty-v ^{-ty)) , a ye E\Eg 

^ 1 y, iiy&Eg 


(A.5) 


This function is an inverse of (p as can be seen be reverting the flow induced by x when needed. 
Also, with the same arguments as before, we can prove that it is C^. 

This implies that cp is a. diffeomorphism from K"* to E. 

Besides, the functions z !->■ Z{z,t) for t > 0, z !->■ and v being C®, p is C® at any interior 
point of {Eg)^. By continuity of for r < s, it can be verified that p is also (7® on the boundary 
dEg. So, <() is a C'®-diffeomorphism from K™ to E. 

Finally, we note that, for any point z^ in dEg^ there exists a point z in dE satisfying : 


\Zg - z = 


< e sup|x(C)l ■ 

C 


xi.Z{z,s))ds 

And conversely, for any z in dE, there exist Zg in dEg satisfying : 


Zg-Z = 


x(Z(z,s))ds 


< s sup|x(C)l ■ 

C 


It follows that 


dH{dEg,dE) < esup|x(C)l (A.6) 

C 

and £ may be chosen as small as needed. □ 

Appendix B. Diffeomorphism extension from a ball. 

Lemma B.l (3.5). Let R be a strictly positive real number and consider a Q diffeomorphism 
p : BrH) -e p{Bii{0)) C K™. For any strictly positive real number e in ]0,1[, there exists a 
diffeomorphism pe ■ 1^"* -7 K"* such that pe{z) = p(z) for all z in cl(B r 10)). 

1 + e 

Proof. It sufficient to prove that our assumptions imply [H Theorem 8.1.4] applies. 
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We let 


U = 

i+f 


A = cl{B^m , 



and, without loss of generality we may assume that ^(0) = 0. 

Then, consider the function F : C7 x / —>■ defined as 

= (1^(0)) ^^,Vte/\{0}, (^(z,0)=z. 


We start by showing that F is an isotopy of U. 

• For any t in /, the function z !->• Ft{z) = F{z,t) is an embedding from U onto Ft{U) C ffi™. 
Indeed, for any pair {za,Zh) in satisfying F{za,t) = F{zb,t)^ we obtain ip{zat) = ipizbt) 
where Zbt) is in The function 'll) being injective on this set, we have Za = Zb which 
establishes Ft is injective. Moreover, we have: 


dz 


(z) 




Vt G/\{0} 


dFo 

dz 


(z) = Id. 


Hence, Ft is full rank on U and therefore an embedding. 

• For all z in U, the function 1 1 —>■ F{z,t) is C^. This follows directly from the fact that, the 
function ijj being C^, and tf){0) = 0, we have 


^{zt) 

t 


dz 


(0)z + z' 



Z^ +o(t) . 


In particular, we obtain ^{z,t) = p{z,t) where 


P(.zA) = ^ 


^{zt)zt - ipizt) 


VtG/\{0} , p(z,0) = -z' 


I , / d’^if) 


2 \ dzdz 


(0) z . 


Moreover, for all t in I, the function z i—>■ ^{z,t) is locally Lipschitz and therefore gives rise 
to an ordinary differential equation with unique solutions. 

Also the set [J(^z t)euxiii^(z,t),t)} is open. This follows from Brouwer’s Invariance theorem 
since the function (z, t) i-G (F(z, t),t) is a diffeomorphism on the open set Uxl. With [T31 Theorem 
8.1.4], we know there exists a diffeotopy Fe from K”’' x I onto which satisfies F^ = F on Ax [0,1]. 
Thus, the diffeomorphism i/je = Te(-, 1) defined on K"* onto verifies ipe{z) = F^{z, I) = F{z, 1) = 

'0(z) for all z G A. D 
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